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ABSTRACT 


Echo splitting is a phenomenon that produces difficulties in 
estimating target range and speed for radar and sonar systems. It occurs 
when more than one returned pulse is received from a single target for 
each transmitted'pulse.' This thesis investigates the effects'of echo¬ 
splitting distortion of a linear frequency modulated (chirp) signal. An 
ambiguity function is proposed for both narrow-band and wideband analysis 
of multiple and overlapping target returns. The double-echo situation is 
examined in detail and equations are developed for both'the'narrow-band 
and wideband cases. Finally, using typical sonar parameters, sample 
computer-drawn plots of these ambiguity functions are presented. 
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CHAPTER 1 


INTRODUCTION 

The recognition of the science of signal processing has led to 
the initiation and development of a class of sophisticated'pulse-type 
waveforms for communication systems. The fundamental property of this 
class is that each member signal possesses a large time-bandwidth 
product. The reception of such a signal by a cross-correlation detector 
produces time compression of the transmitted pulse. Selection of 
suitable transmitted waveform structures is called "signal waveform 
design." Systems utilizing correlation detection for signals having 
large time-bandwidth products are commonly said to employ "pulse 
compression techniques." 

A pulse-compression technique that has-found; wide acceptance 
in radar and sonar applications is based upon a. linear frequency 
modulated (LpM) pulse. In this technique, the carrier frequency is 
linearly swept to form a transmitted pulse of constant'amplitude. The 
instantaneous frequency is expressed then as 

fo * 

during the pulse,where f^ is the initial or carrier frequency. See 
Figure 1-1(a). 

One method for implementing cross-correlation detection at 
the receiver is to use a linear filter whose transfer response to a 
unit impulsive input is the time inverse of the transmitted signal 
pulse. This is called a "matched filter detector," For LFM, the 
matched filter should have a constant gain characteristic over the 
required passband and a phase characteristic designed to provide a 
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linear time delay versus frequency plot. Figure l-l(b). By retarding 
the leading portion of the received pulse waveform more than the trailing 
portion, the net effects of matched filter detection will be a compres¬ 
sion in the time dimension of the filter output and an increase in the 
received signal amplitude.^ 




(a) Transmitted Pulse (b) Matched Filter 

Fig. 1-1 

LFM is both conceptually and technologically simple. The efficiency of 
operation of LFM systems has resulted in the continuing employment of 
this technique in advanced radar and sonar applications. 

2 

The ambiguity function, originally defined by P. M. Woodward, 
examines the relative merits of a particular type of signal processing. 

It is a measure of the ability of a signal processing scheme to resolve 
two targets differing in range by a time delay t and in doppler shift by 
frequency f'. But the ambiguity function can be interpreted in another 
manner. Rather than investigating the resolution between two targets, 
a replica of the transmitted signal is used at the receiver to simulate 
one target. Then the received signal is cross-correlated with this 
locally produced replica signal. W. M. Siebert has indicated that this 
cross-correlation of the received signal with a locally produced replica 
signal can be considered as the output of a matched filter. Consequently 
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the ambiguity function can be thought of as the envelope of the output 
of a matched filter, in the absence of noise, for a signal differing in 
range delay by a time t and in frequency by f. 

It is the intent of this paper to investigate.the ambiguity 
function of a linear FM signal in the presence of echo splitting. Echo 
splitting is a broad classification of two independent phenomena that 
cause distortion in the received signal. 

The first of these is attributed to multipath transmission 
between source and target. Multipath manifests itself particularly in 
sonar systems where bottom bounce, surface reflection, and direct path 
propagation, and various combinations of these paths can illuminate the 
target. Figure 1-2 illustrates these paths by which illumination of 
the target can occur. 



Multipath is also found in certain radar applications and is especially 

4 

noticeable in airborne radar. In this instance, a target can be il¬ 
luminated by direct path and by surface reflection as shown in Figure 
1-3. 
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source 



The vagaries of both surface, and bottom reflection.wtll provide many 
paths between source and target. For this study, discussion will be 
limited to four paths: 

1. A signal sent and received over a direct path. 

2. A signal sent and received over a'reflected path. 

3. The path consisting of transmission over the direct 
path and reception via the reflected path. 

4. The'path consisting of transmission over'the reflected 
path and reception via the direct path. 

In the sonar case, the time delay difference between a signal 
sent and received over a direct path, .a signal sent and received via 
bottom bounce and a'signal that takes one path to’the~target and the 
other for its return to the source may often* be sufficient:to resolve 

5 

the returned echos. However, the signal that Is*sent’direct path and 
returns via bottom bounce will arrive In close proximity to a signal sent 
to the target via'bottom bounce, and returned direct path. In a sonar 
system, separation between these two paths on the order of a few tens 
of milliseconds can be expected. Time separations of this order will 
produce difficulties in the detection problem. 

The second cause of echo splitting is peculiar to'sonar systems 
and is due solely to target dimensions. A sonar target is normally of 
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sufficient length, when compared to the length of a pulse, to display 
multiple reflecting surfaces. These individual reflection points will, 
depending on the aspect, present returns to the'receiver sufficiently 
separated in time so as to cause distortion. The detection problem is 
further complicated by the possibility that'the'reflecting points can 
display different doppler shifts as in the case'Of:a turning target. 

There have been many developments of both narrow-band and, 
more recently, wideband ambiguity functions.' In the narrow-band case, 
the assumption is made'that the only effect of doppler on a*transmitted 
signal is to shift the frequency components, of the signal. This assump¬ 
tion neglects the effects of doppler distortion on the envelope or 
modulation. These effects will be shown to have considerable importance 
in wideband analysis. 

Chapter II contains the development of a general narrow-band 
ambiguity function for the echo-splitting case. A'model'that explains 
how this function could be calculated is proposed. Finally, a specific 
example employing two returning signals is derived and computer-drawn 
graphs of this function are presented using typical'sonar parameters. 

In Chapter HI, wideband'considerations are introduced, and a 
simplified derivation of these effects on'the ambiguity function is 
given. 

Chapter IV proposes a system model based upon the wideband 
ambiguity function. An example with two returning signals is presented 
for the wideband case, and computer outputs'are shown for various delays 
between the two signals. 

Finally, Chapter V'surtmarizes the results of this'investigation, 
and suggests possible extensions of the work conducted in this thesis. 
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CHAPTER II 


NARROW-BAND ANALYSIS 


A. BACKGROUND 

As previously mentioned, the work of Siebert best suits the 
initial phase of this investigation. Consequently, the following 
development is largely taken from that paper and is intended to clarify 
the symbols and expressions to be used. 

In a linear FM system, the transmitted signal, s^Ct), can be 
written as 


So(t) = A{u(t) - u(t-T)} cosCco^t + j yt^). 


(2-1) 


where 



u(t) = the unit step function = 


For ease of expression, we define 


X(t) = u(t) - u(t-T), 
so that Equation 2-1 is written 


s^(t) = A X(t) cosCwpt + i yt^). 


It is customary to normalize the energy level of the outgoing 
pulse by stating the condition that 



( 2 - 2 ) 


Carrying out this operation reveals that the amplitude A of the pulse 



should be chosen 


The transmitted signal can also be written ina canonic form, 
in this case as the real part of a complex signal. Thus 


So(t) = Ris^(t) 


where 


syt) . X(t)|S^(t)| 
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The complex modulation factor S^{t) must conform to the normalization 
of Equation 2-2: 

|So(t)| = A . 

In the narrow-band case, a returning signal, delayed by an amount t and 
frequency shifted by can be written as 

s (t-T) = Re[S(^ (t-x) , 

'*^d d 

where 


S (t-T) = a X(t-T) |S . 

“d 

In this last expression, a is introduced to express the amplitude of 
the received signal as a percentage of the normalized transmitted signal. 

In the case of echo splitting, the received signal is a sum of 
the N returning echoes. It can be written as 


r(t) ^ il a-s (t-t : . (2-3) 

j=l dj j 

where t. is the range delay, w. is the frequency change, and a. is 

j A 

the relative amplitude of the reflector. 

It is necessary at this point to digress somewhat in order to 
give a short justification for the present approach. At'first glance, 
it appears that the normal method'of testing to see if'two returning 
signals can be resolved wi11 prove not only simpler'but also more 
enlightening. For example, in the case of a' target with two returns, 
the usual approach is to compare' the ambiguity function of one return 
with that of the other return. If the two curves intersect below the 
3db points, then the two targets are said to.be resolvable. That is, 
an operator should be able to distinguish two targets. 

But this is not the problem'treated'in this paper. The multi¬ 
ple returns that arrive with different delays, hence different phases, 
add together in the matched filter. Therefore, if an ambiguity function 
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can bfindeveI oped utilizing these considerations, it could possibly sug¬ 
gest means of target identification and recognition. Such a possibility 
argues strongly for an effort in this direction, particularly in the 
sonar environment. 


B. DEVELOPMENT OF THE NARROW-BAND AMBIGUITY FUNCTION 


Siebert developed- the- relationship' between the- a-posteriori 


probability of detection, P(A, x, to), and the a-priori probability. 


Pq(A, t, (o), as 


P(A, T, to) = K Pq(A, t, to) e^o 




“i e "o (2-4) 

For a single echo, x and to represent the range delay and doppler fre¬ 
quency mismatch between the-received signal r(t) and "the replica of 
the transmitted signal s'^ (t-x). N^ is the noise power density in 
watts/Hz and K is a normalization constant-. It'is obvious from 
Equation 2-4, that’ the integral 


A “ I (t-x) dt (2-5) 

contains all the information that is necessary in'the'detection'process. 

It is also obvious that with the' proper selection of the factor 
s' (t-x). Equation 2-5 is simply the output of a matched'filter. This 

(Oi 

integral will, however, contain the' rapid variations of the fine struc¬ 
ture' of the signal caused by the carrier phase term, w^x. The fine 
structure does not in any way aid in the determination of'the range of 
the target. For, although-it is possible to*determine'in which fraction of 
the wavelength the target Is located, there is no. way to predict at 
which particU'lar half-wavelength this occurs. 

Siebert solved the problem in the narrow-band case by treating 
the (jj^x term 1n s'^ (t-x) as a random variable, 9, uniformly distributed 
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in the interval 0 £ 0 £ 2 tt. Performing this operation on Equation 
2-4, he arrived at the expression 

A" 

- Tjm— 

iP(A, T, oi) * KP(A, T, w)e o * 

o 

where 

Iq is the modified Bessel function of the first kind and zero order 
and 


A = 


/ ^(t)S^ (t-T)e'^“o^ dt 


(2-6) 


This operation has the result of essentially taking the envelope in the 
T domain, of the expression In Equation 2-5. Equation 2-6 will be used 
as the basis for defining the ambiguity function ip. For a multiple echo 
return, let 


ip(T,a)) = 


I 


r(t)S'^(t-T)e 




dt 


where now t and w are multidimensional and must be carefully interpreted. 

Equation 2-6 differs fronr Siebert's ambiguity function in that 
the carrier frequency phasor, -fs retained. The necessity for 

keeping the carrier frequency can best be explained by consideration of 
a continuous-wave (cw) pulsed transmitter system. If two point reflec¬ 
tors are separated, for example, by a distance then'the received 
signal would essentially vanish since the total "path difference 1s 
Similarly, for a separation of resulting In a total path difference 
of X, the received signal would be doubled. Consequently, we might 


17 










expect that 


ii<(0,0) = 0 for a total path difference of ^ 

i|;(0,0) S 2 for a total path difference of X. 

It is possible, to hypothesize an equation' to'meet' these conditions. 
By defining the time' for the'wave to'travel?the'roundftrip distance 
between the two reflectors as.A, then for 
cA - 



the result is 


ij;(0,0) = 0. 

Also, for 

CA = X 

A - 

A = 77" • 

then 

ii;(0,0) = 2. 

It is easy to see that 


t |»( 0 , 0 j '='^1 


UqA 

cos -|- 


satisfies these conditions. Appendix' Ai shows that for the case of two 
point reflectors'," Equation’ 2'-6 .1 eads? to- this'expression: for i|j(0,0). 

The carrier frequency phasor 1 s-' thus "an' important factor'! n ambi gui ty 


functions-pertaining’ to' signalsrin'ar echorsplittingrenvironment. 

Before continuing',“itris"appropriate tn'indicate just how this 
function- might be' generated in"practice-, fpigure 2-T is? one’ possible 
block diagram- for a rratched? filter system that is- suitable’ for narrow- 
band considerations’. Figure? 2-2 is a breakdown of a'particular doppler 


18 










RepUca 

Signal 


si(t) 


Frequency 

Shifter 


s' (t) 

Wi ' ' 


Received 

Signal 


r(t) 


(t) 


(t) 

n 


Doppler 

Filter 

Banks 


> V'(x,w,) 

V 


} 

O 

o 

0 


? 

o 

« 

—^ 



Display^ 


Fig. 2-1 ENVELOPE OF THE OUTPUT OF A MATCHED FILTER 
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banko The functional blocks entitled "collator" collect the various 
outputs for each shifted value of t and w and arrange the values in 
order on the x and w axes' respectivelyo 

C„ DOUBLE-ECHO EXAMPLE 

If only two reflectors are considered^ then the received 
signal, T(t), can be written as 


r(t) - a, s (t-t^ ) + cx.^s (t-t^^)„ 


'd, % 


(2-7) 


Substituting this expression Into Equation 2-6, 




r 



By defining 


t' = t-t 



and t“ = T-t 


d, ’ 


Equation 2-8 can be written as 



Since 


gjwotd, _ ^ ^ 


this becomes 





'd 


(t) + a.s (t-A) S' (t-T')e'^‘^o^ dt .(2-9) 

^ (Jl)j I tjO» 
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In Equation 2-9, the dummy variable of integration t' is replaced by 
t. From the above substitutions, it is seen that A is simply the 
difference in arrival time between the two returns. Thus, Equation 
2-7 could have been written as 

>■(») “ “.V 

d, dg 


although it should be noted that the ambiguity function is now plotted 
as a function of t' which is shifted from the original t origin by the 
time delay of the first return. 

For a linear FM signal. Equation 2-9 is written as 


ij)(T',a),A) = ^ + OgXCt-A) 

.cos{(a,„ + o.j^)(t-A) + 


(2-10) 


By applying appropriate trigonometric identities, Equation 2-10 can be 
written in the form 




“t ^ 

,a),A) = [ ^0 


cos X + cos y + jsi 


n X + j sin^^ 


dt 




cos z + cos V + j sin z + j 


j sir^^ 


dt 


. ( 2 - 11 ) 


where 


X = yt^ + (2a)Q + + Wj - yT')t + %yT'^ - 

y = (o) - oj. - yT')t + - 0 ) t' 

X u J X 

z = yt^ + {2a)^+a)j^+a)j - y{T'-A})t + %y(T'^+A^)-(a)^+a)j )A - 
V = “ y{T'-A})t + )A + %y(T*-A^) + . 
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The limits of the integrals are determined by the envelope factors 
and the particular value of t'. 

These expressions can be significantly reduced by employing 
narrow-band assumptions. Integrals of the form 


Jcos(ax^+2bx+c) dx 'J^^cos 


C 


+sin(^) S 

a 


[p 


(2-12a) 


and 


Jsin(ax^+2bx+c) dx = 


cos(^) SI 


-sln(^) C 

a 


Ijil' 

P 


(2-12b) 


(ax+b) 


can be shown to be negligible under narrow-band considerations. 

In these expressions, C( 2 ) and $( 2 ). are the. Fresnel integrals 
2 


0 ( 2 ) - J cos(^t^) dt 

O 

Z 

S(z) = j sinC^t"^) dt. 


6 


Introduce the argument: 
„ b^-ac 


(2-13) 


If only the first reflection is considered, then 

a = p 

b = % (2w +w,+co. -yx') 

o 1 u, 

C - 


% yx' 


For 03^ and « o)^. Equation 2-13 reduces to 


a3‘-ya3 
o o 


(2-14) 
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If typical values of carrier frequency and rate of radian fre¬ 
quency sweep are substituted in Equation 2-14, it is not difficult to 
see that the expression Is relatively independent of t'. Hence, we 
may write 

0) 2 

0 S —2— • 

P 

The next item to be examined is the argument of the Fresnel 
integrals 



Again under narrow-band considerations', this expression can'be written 
approximately as 



Not only is this term relatively independent of t', but also the argu¬ 
ment is a large number for representative values of carrier frequency 
and rate'of radian frequency sweep. The values of both the sine and 
the cosine Fresnel integrals approach % for large arguments. Conse¬ 
quently, the right'hand member of Equations 2-12a and 2-12b can be 
written respectively as 

_ I 

cos e + sin e) 

and ' 


^^(cos e - sin e) 


For a - p, these expressions'contribute'little to the integral in 
Equation 2-11, Therefore', terms having the form of Equations 2-12a 
and 2-12b can be neglected. 
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Hence, Equation 2-11 is written as 


ij;(T' ,w,A) = 




yi ')t + k yx' 


A 


Wd -y{T'-A})t + (wo+cod )A + % y(T'^-A^)-WiT'] 


dt 


which can be reduced to 


i|;(t' ,co,a) = 


a 

_lJ gdCw^-wd^-yx 

T n 


-yx')t 


dt 


(2-15) 


Ah 

+ J gj[(wi-wd2-y{'r'-A})t +(a)Q+wd2)A - % yA^] 

3 


Equation 2-15 was coded in Fortran and programmed on an IBM 360 
computer. By assuming a sequence of equispaced values for target doppler, 
a family of curves of 4; versus x were computed for the zero doppler filter 
bank. The output was plotted by the computer in the form of an oblique 
projection as shown in Figure 2-3. The origin of the graph scale on the 
X axis is set in relation to the particular curve at which doppler match 
occurs. In this figure, the middle curve represents doppler match 
between filter bank and target. 

In order to affect'comparison with'the wideband ambiguity func¬ 
tion derived and plotted' in Chapter IV, sonar parameters: were assumed. 
Consequently, the long pulse lengths common to'sonar systems restrict 
the amount of pulse compression that can be assumed if the signal is to 
be considered narrow-band. (The requirements for a signal to be con¬ 
sidered narrow-band are given in the next chapter.) 
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The parameters of Figure 2-3 are: 

f = 1000 Hz 
0 

k = 10 Hz/sec 

T = 1 sec 

A = .004 sec 

with equal echo amplitudes 

a = a =1. 

1 2 

Figure 2-4 is a plot of the same function for a delay of' .OT seconds. 

It is interesting to contrast' Figure 2-3, the result of applying 
Equation 2-15 to a narrow-band signal, with the results'obtained when 
a wideband signal is assumed and substituted in to this same narrow- 
band equation, f 2-5 is a plot for this wideband signal, which 

has the same parameters as the signal used in Figure 2-3 except that 
k = 200 Hz/sec. 

In this graph, the delay A is distinctly observable*. The echoes almost 
cancel, because of the wide bandwidth; Also*noteworthy"is the relative 
invariance of the doppler cuts of the ambiguity function over the com¬ 
plete range of target doppler. 
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CHAPTER 3 


WIDEBAND AMBIGUITY FUNCTION 


A. BACKGROUND 

As previously stated, the implications inherent 1n the use of the 
narrow-band'ambiguity function are that the doppler effect does not 
significantly alter the envelope of the transmitted'signal. This is a 
valid assumption in most radar cases. However, for sonar systems, the 
combination of lower carrier frequency, longer pulse length, and slower 
medium propagation speed makes this assumption restrictive. 

The commonly accepted cross-over point between the narrow-band 
case and the wideband case occurs when^ 


where 

c = propagation speed in the medium 
B = signal bandwidth 
T = pulse width 

v^= radial velocity (assumed positive when the range is closing). 
It has been shown that for the case of non-accelerating targets 
the inequality can be confined to^ 



Taking representative sonar figures of 
B = 200 Hz 

T = 1 sec , 

we find that 

''r ~ 1.5 knots. 

It is obvious that narrow-band considerations do not apply. 
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Wideband ambiguity functions have been under Investigation for 
some time. Perhaps the most noteworthy treatise on the subject is a 

g 

paper by E. 0. Kelly and R, R„ Wishner, which presents a very detailed 
mathematical development of a wideband ambiguity function. References 8, 
9, and 10 investigate the effects of an accelerating target on linear FM 
systems. Finally, References 11 and 12 expand'the concept of a gener¬ 
alized ambiguity function. 

S. A. Kramer utilized frequency versus time plots to examine the 

9 

effect of wideband analysis on linear FM. His method offers the best 
insight into the problems associated with"wideband systems. Therefore, 
the following section is a brief synopsis of Kramer's development. 

In the narrow-band case, when there is a doppler mismatch, the 
received signal is simply shifted in frequency uniformly throughout its 
length some distance from the replica signal.: The lines remain parallel 
because the pulse compression factor of the'received'signal is the same 
as that of the transmitted signal. 



Fig. 3-1 Mismatched signals: Narrow-band case 

In the wideband case, however, a slope change is. introduced, 
caused by a pulse compression factor different from the pulse compres¬ 
sion factor of the transmitted signal. Thus, if the replica signal does 
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not match the received signal in doppler, the two lines are not parallel. 


f 



received 

signal 


1-t 

Fig. 3-2 Mismatched signals: Wideband case 


Of course, if the received signal is.exactly'matched in doppler, 
then'the'two lines would be parallel even for the wideband case. 

If the discussion is limited to targets of ^constant velocity, 
there are'three effects of doppler on'the returning signal: 

1. The carrier frequency f^ of the transmitted signal is 
shifted by the doppler frequency. 

2. The envelope of the transmitted signal is altered. 

3. The energy of the received signal differs.from the trans- 


13 


mitted signal energy. 


In the narrow-band case only the first effect is taken into consideration. 

The following analysis reexamines the frequency shifting effect 
and introduces the other two considerations into the development of a 
wideband ambiguity function. 

Each frequency component of the received signal, can be written 


as 



(3-3) 
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where 


= received frequency component 
= transmitted frequency components 
For v^« c, Equation 3-3 can be approximated as 

f, :(i+^)ft. (3-4) 

It is convenient to define a doppler factor, D, where 
2v 

D = 1 + 

c 

The doppler frequency, f^, can be defined as the change in the 
initial frequency of the transmitted signal, fo^» or 
fd' fo,- for 
But, from Equation 3-4, 
fOy-“ 

fj = (D-1) fot, 
and 



0 


which shows the relationship between the doppler factor and the doppler 
frequency. 

The transmitted signal can be written in the frequency domain by 
taking the Fourier Transform of s(t), 

S(f) = F[s(t)]. 

The received signal can be expressed in terms of the transmitted signal 
in the frequency domain as 
R(f) = 
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Now the received signal can be written as the inverse Fourier 
transform of S{^), 


r(t) = F-‘[S(^)] = / S(^ df . 

—00 

Introducing a change of variable, f ^ » Equation 3-6 becomes 
r(t) = S(f')e'^^'^^'^^'‘^^ df . 


(3-6) 


(3-7) 


This is the normal inverse Fourier transform, except that the trans¬ 
formation is into the revised time domain Dt. Equation 3-7 becomes 
r(t) = Ds(Dt), 

which for a linear FM signal is written 

r(t) = ADX(Dt) cos + % y(Dt)2}. (3-8) 

If Equation 3-8 is now examined, it can be seen that the three 
effects of wideband analysis have occurred. 

1. The frequency components of the transmitted signal are 
shifted by the doppler frequency. For, the. instantaneous frequency of 
the received signal can be written as 

uJqD + pD^t . 

But from Equation 4-5 this can be written 


(Do(l + :p) + pD^t 

0 

= (*) + u). + pD^t 
o d 

Thus, not only has the initial frequency been shifted but effectively 
there is a new pulse compression factor pD^T. This accounts for the 
difference in slopes of Figure 3-2. 

2. The envelope of the transmitted signal is altered. This has 
occurred through the new definition of the envelope term X(Dt). 


33 




For 


X(Dt) ^ u(Dt) - u(Dt - T) 

= u(D{t}) - u(D{t - J}) 


- u(t) - u(t - . 


From this expression., the expansion or compression of'the pulse is 
shown to be dependent on the value of:the doppler factor D„ 

3, Finally, the energy of the return signal has been changed 
by the doppler factor D;, Reference 13. showed; that' the* ratio of the 
received signal energy density to the transmitted signal energy density 
is given by D„ 

The energy of the signal'of Equation 3-8 


J°A2D2X(Dt)X(Dt) 


cos^(w Dt:+ \ pD^t^) dt 
o 


T/D 



o 


" 2 J 0 

= D. 


B. GENERAL DEVELOPMENT OF THE WIDEBAND AMBIGUITY FUNCTION FOR ECHO 
SPLITTING 

The output of a matched filter is written as 



(3-9) 


which will be recognized from the narrow-band discussiono 

Investigation of Equation 3-9 might be accomplished by writing 
the received and replica signals in canonic form as: 
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r(t) 


= Re[R(t)e^“o^] 


(3-lOa) 


s'a)j^(t-T) = Re[S'a)^(t-T)e^‘*^o^^"^^] , 


where as before. 


R(t) 


- DAX(Dt)ei«“<f* ^ ^ 


and 


S'w^(t) = X(Dt) 


S'o(t) 


JKt + ig uDH^) 


Recalling that Re[ 2 (t)] can be written as 
%[2(t) + 2*(t)] , 

Equation 3-9 can be written as 


(3-lOb) 


^ = h J°[R(t)e^*“o*S'a)^(t-T)e^*“o^*"'^U R(t)eH*S w^(t-T)e‘H^*’'^^ 

Joo 

(3-11) 

+ R*(t)e‘'^*“o*S'Wi(t-T)e'^'“o(''^"'^UR*(t)e"'^“o*s'i^(t-T)e"^“o^*“'^^] dt. 

If the terms, R(t)e^“o^ and S'u)i(t-T)e'^“o^^“'’-^ are analytic, 

then since it is possible to state that the correlation integral of two 

analytical signals or the correlation integral of the complex conjugates 

11 

of two analytical signals iszero. Equation 3-11 reduces to 

A - Jg j*R(t) S'wj^(t-T)e’^*“o'^dt. (3-12) 


As this is the output of a matched filter, the ambiguity function can 
be arrived at by simply taking the envelope of Equation 3-12. If such 
an envelope could realistically; be formed (in the wideband case, this is 
not a moot question), * the ambiguity function would be written as 


i|^(T,a)) 



R(t) s'S^(t-T) dt 


(3-13) 
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The fine structure, attributable to the carrier phase term e'^“o^, has 

been omitted.by the process of taking the absolute value. It is the 

same function as developed by Siebert in the narrow-band case. For 
narrow-band analysis, this thesis rejected Siebert's'definition of 
the’ambiguity function in the echo-splitting environment, because of 
the absence of the carrier frequency term. In addition, the assump¬ 
tion that' an analytic signal model is appropriate under wideband condi¬ 
tions is highly suspect. 

Consequently, the only recourse is to accept Equation 3-9 as the 
output of the matched"filter and begin at that’point to'devise a suitable 
ambiguity function. As it is the envelope of this expression that is 
being used for the definition of the ambiguity function, the fine struc¬ 
ture is removed simply by taking the absolute value of Equation 3-9. 

The ambiguity function is then written 





Wj^(t-T)e'^''^o^dt 


(3-14) 


which is the same form as that used in Chapter 11 for the narrow-band 
case. 
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CHAPTER 4 


WIDEBAND AMBIGUITY FUNCTION FOR DOUBLE-ECHO CASE 
A. SYSTEM MODEL 

If the time origin is taken as the starting point of the 
transmitted signal, then Equation 3-8 may be written as 
r(t) = ADX(D{t-t^}) cos{a)oD(t-tj) + % yD2(t-t^)2} . 

The system model for generating a wideband ambiguity function might be 
constructed to be similar to Figure 2-1. Figure 4-1 and Figure 4-2 
are'possible models for the wideband analysis. 

The output of the matched filter for the range estimate 
T|^, and the i^^ doppler filter, is written 

y(Tk,Wj.) - Y S cos{(o^D'(t-t^j) + % ]iD^(t-tj)^} 

• X(D^t-T|^) cos{w^(D^t-T|^) + % p(D^t-T|^)^}dt. (4-1) 

Defining 

f = t-tj 

and 

Equation 4-1 is written as 

y(T',w ) = I J DX{Dt') cos{w Dt' + %i(Dt’)nX{D (t'-T')} 

* .00 ® ^ 

• cos{a)^(D^[t'-T'])+ %ij(D^[t'-T'])^>dt'. (4-2) 

This can now be referred to the time base of the replica signal by 
defining 

•k 

D^t' = t 

and 

D = D'. 
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Fig. 4-1 ENVELOPE OF THE OUTPUT OF A MATCHED FILTER 




y(T,w^) 



To 0 ) collator 


Fig. 4-2 DOPPLER FILTER 
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With these substitutions. Equation 4-2 can be written as 


y(T',w.) = f J D'X(D't*) cos{a)^D't* + %y(D't*)2}X(t*-D^T') 

^ ' -00 ^ 

“ cos{w (t*-D,T') + %y(t*-D.T')^} dt . (4-3) 

0 1 1 

From the substitutions that have been made, it is seen that when 

D.t' = 0, 

1 


Thus, if the system of Figure 4-1 is used, the doppler factor enters 

into the determination of range. 

If, on the other hand, another model for the system is chosen, 
this problem does not arise. Figure 4-3 and Figure 4-4 describe such 
a system. For, with this as the model, the output of the matched filter 
for the range gate and the i^^ doppler filter is written as 




( DX(D{t-tj})cos{w^D(t-tj|) + 


• X(D.{t-T,}) 
1 k 


cos{a)^D^(t-T^) + Jgy(D^[t-T^])^}dt. 


(4-4) 


Now if we define 

t' = t-t^ 
and 

" V^d’ 

Equation 4-4 becomes 

y(T',o),) = Y ^ DX(Dt') cos{u^Dt' + %y(Dt')^}X(D^[t'-T']) 

-00 

• cos{a) D.(f-T') + %y(D.[t'-T'])^}dt'. (4-5) 

O 1 1 

As in the preceding example, this expression is now referred to the time 
frame of the replica signal and is written as 
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Fig, 4-3 ENVELOPE OF MATCHED. FILTER OUTPUT 
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Fig. 4-4 DOPPLER BANK 
































y(T'.a) ) = fo'XCD’t*) cos{w D't* + %y(D't*)2}X(t*-D x') 

I Joo ^ ^ 

• cos{w^(t -D^t') + JgyCt -D^t')^} dt*. (4-6) 

Now when 

D^t' = 0 , 

it is seen that 

T, = t,. 

k d 

It is for this reason that Figure 4-3 is chosen as the model of the 
system. Accordingly, Equation; 4-4 describes the output'of the matched 
filter. 

B. SINGLE-RETURN DEVELOPMENT 

It is interesting'to compare'the'ambiguity function that was 
derived in this paper with those'previously developed.. As mentioned 
earlier, the essential difference in our proposed ainbigutty function 
is the inclusion of the carrier frequency term This term was 

included because of the assumed importance of the carrier frequency on 
the ambiguity function in the echo-splitting case. 

^ The following section develops an expression: for the'ambiguity 
function arrived at in Equation 3-14 in the case of linear FM when only, 
a single return is experienced. 

Equation 4-6 can be written as 

y(T',w^) = Re| 

As in the narrow-band case, we have defined the"ambiguity function as 
the envelope of the output of a matched filter and we approximate this 
envelope by taking'the absolute value: of the expression in'Equation 4-7. 


I D'X(D't) cos{w^D't-+ %ii(D't)2}X(t-D^T') 

gj[(Oo(t-DiT') + %y(t-DiT-)^] 


(4-7) 
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The final result for the linear Ff4 case with single return is 



^{t',u)^) = 7 J D'X(D't) cos{u)^D't + %y(D't)2}X(t-D^T') 

I _oo 


(4-8) 


Computer outputs of this function utilizing various sonar 

parameters are shown in Figures 4-5 through 4-7. One interesting 

property of this function is the apparent splitting of the pulse for 

the non-matched condition. It is also interesting to note the narrow 

doppler tolerance of this function. This is true of all wideband 
g 

functions. 

C. DOUBLE-ECHO EXAMPLE 

If the echo-splitting case is now considered, the received 
signal, r(t), is written as 


N 



Assuming once again that there are only two reflectors. Equation 4-9 
becomes 


r(t) = Aa,D,X(D,[t-tj_]) cosCw^D, (t-t^ ) + [t-t^ 


+ Aa2D2X(D2[t-td23) cos{w^D2(t-td2) + ^pCD.Ct-tdg])^} . 
y(T,u^) is then written as 
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FIG. 4-5 SINGLE-RETURN WIDEBRND RMBIGUITV FUNCTION 
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FIG. 4-6 SINGLE-RETURN WIDEBRND RHBIGUITV FUNCTION 
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If the envelope of this expression is now taken. Equation 4-10 is 
wri tten 

=y| cos{a)gD,(t-td^)+%y(D,[t-td ])n 

+ a2D2X(D2[t-td^]) cos{woD2(t-td2)+Jju(D2[t-td2])2}J (4-11) 

By defining 

f = t-td^ 

A = 

and 

t' = T-td^ , 

Equation 4-11 becomes 






a,D,X(D,t')cos{wQD,t'+:!gu(D,t‘)^}+a2D2X(D2[t'-A]) 


• cos{WoD2(t'-A)+%ii(D2[t'-A])^}J (4-12) 


Finally, by defining the quantities 



and referencing all returns to the replica signal by the transformation 
D^t' = t* = t , 


47 






Equation 4-'i2 can be expressed as 



a , D' X (D ■ t) cos {(jOqD ' t+%y (,D ‘ t) ^ }+a^ D"X (D" [t-D^A]) 


•cos{a)oD"(t-DiA)+%vi(D"[t-DiA])2}] 

‘X(t-DiT')eJ*^“ot+%M(t-DiT')^] 


(4-13) 


As a final substitution we define 


The final expression for the ambiguity function of the linear FM system 
is 



+ a^D"UD'' [t-D^A]) cos {u^D" (t-D^A)+%y (D" [t-D. A]) 


(4-14) 


It is important to note a fundamental difference between the 
ambiguity function above and that developed in the narrow-band case. 

In the narrow-band case two paths are open for further investiga¬ 
tion: 

1. Plots of different replica signals can be correlated with a 
target of particular range and doppler. 

2,. Plots of a particular replica signal can be correlated with 
targets of various velocities. 

In the wideband case, the fact that our function is plotted versus 
t", which by definition is proportional to D^, restricts the choice to 
Case 2, if meaningful results are to be obtained. 

It can be shown that with suitable approximations, Equation 4-14 
reduces to Equation 2-10 developed in the narrow-band case. Appendix B 
develops the various equations needed to evaluate Equation 4-14. 
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Figures 4-8 and 4-9 are sample plots of Equation 4-14 for the 
0 doppler bank. The plotting procedure is the same as that used in the 
narrow-band analysis. Basic parameters for these plots are: 

f = 1000 Hz 
0 

k = 200 Hz/sec 

T = 1 sec 

and 

a, = ttj = 1. 

For Figure 4-8, a delay of .003 seconds is assumed and for Figure 4-9, 
a delay of .004. No significant difference exists between these two 
families of curves. 

Figure 4-7 may be compared with Figures 4-8 and 4-9. Since 
Figure 4-7 represents results for a single echo, = 1 and * 0, 
one must mentalTy double its ordinates before comparison. Furthermore, 
the doppler intervals between curves are 2.5 times greater. With these 
qualifications, the figures show that for A = .003 and .004, the ampli¬ 
tude of the ambiguity function is less than for A = 0, and the peak is 
somewhat broadened. With a time-bandwidth product TBof 200, the main 
peak for the zero doppler cut of the single-echo ambiguity function 
should have first'nulls at t *= + ^ = + .005 sec. The limiting value 
of A for double-echo resolution is thus approximately 5 milliseconds. 
For A = 10 milliseconds, a complete echo split should occur. 

Time limitations, as required for computer plotting, precluded 
generating further echo-splitting exhibits at this time. 
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FIG. 4-8 WIDEBAND GRAPH, DELRV=.003 
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FIG. 4-9 WIDEBRND GRRPH, DELRV=.004 











CHAPTER V 


CONCLUSIONS 

A. SUMMARY 

The purpose of this thesis was to formulate an ambiguity func¬ 
tion of a linear frequency modulated signal under both narrow-band and 

i 

wideband conditions in order to investigate the effects of echo splitting. 

Because echo splitting is primarily a problem in underwater detection 
and'tracking systems, representative sonar parameters were utilized. 

Computer-drawn plots of the ambiguity function for the case where two 
returning signals are received were presented. 

I 

It was suggested that when echo splitting is caused by a target 
exhibiting multiple reflecting surfaces due to the physical extensions 
of the target, there may be a possibility that the'multiple returns 

I 

could provide a means of target recognition. The graphical results ; 

i 

obtained during'this investigation were too limited to evaluate the | 

effects of echo splitting on target signatures. 

B. AREAS OF FUTURE STUDY 

Further investigation of the wideband case appears warranted 
because of its particular importance to sonar systems. It would also 
seem profitable to investigate the effects of echo splitting on other 
large time-bandwidth product type signals such as Barker codes, pseudo¬ 
noise codes and complimentary codes. < 
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APPENDIX A. 


EXAMINATION OF THE NARROW-BAND AMBIGUITY FUNCTION FOR A PULSED CW 
SYSTEM UNDER MATCHED CONDITIONS 


For vi=0, Equation 2-10 can be written as 
i|^(T,a),A) ^ I"! cos{(a)Q+a)(j^ )t} + a^XCt-A) 

cos{(a)Q+a)j|^)(t-A)}j X(t-T')e'^'^‘*^'^ ^ ^o^^ dt 


. (A-1) 


Consequently, ip(0,0,A) can be written, with some trigonometric substitu¬ 
tions as 


11^(0,0,A) = y| j<^cos(2wQt) + 1 + j sin(2a)^t^ 


dt 




COS(2a) t-oj^A) + COSw^A + j sin(2a).t-(ji)^A) (A-2) 


0 0 


0 0 


+ j sinw 


> 


where 


a, = = 1. 


Equation A-2 can then be written as 

i|;(0,0,a) = 1 r2cos^a) t dt + cosu A../"zsin^o) t 

I J ^ ^ A ^ 

[o A 

+ sinw^A rsin2a)Qt dt + j sin2a)_t dt 


dt 


^sin2a)Qt dt + sin2a)Ql 

+ sinw^A ^ 2sin^a)Qt • 


(A-3) 


If the assumption is made that there are an integral number of cycles 
within the envelope of the pulse, 

Y = T; n = 1,2,... , 
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th6n considGfdble rGduction 1n Equation A-3 is possible. For example, 


^ 2cos^i 


aj_t dt = T 
0 


sin 2w^t dt = 0. 


As this investigation is directed at values of A on the order 
of a few milliseconds in the sonar case, the terms 

yT 


cosWqA ^ 2sin^w t dt = -TcoswqA 


and 


cosw^A f sin2a)_t dt = 0. 

0 1 ° 


These same arguments hold for the imaginary terms and so Equation A-3 
can be written as 
i|>(0,0,A) = 

1 + 


1 + cosWqA + j sinw^A 


= 2 


wqA 

COS 


(A-4) 


*The value of this integral will actually be slightly less, due to 
the process of integrating over a partial cycle and reduction by A. 
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APPENDIX B 


EVALUATION OF THE WIDEBAND AMBIGUITY FUNCTION 

Integration of Equation 4-14 is facilitated by investigation of 
the envelope terms. X(t) was previously defined as 
X(t) = u(t) - u(t-T). 

Therefore 

X(D't) = u(D'{t}) - u(D'{t-^,}) 

= u(t) - u(t-^,). 

In a similar manner, 

X(D"{t-D.A}) = u(t-D.A) - u(t-D.A-^„) 

and 

X(t-T") = u(t-T'‘) - U(t-T-T"). 

The received envelopes can be plotted as 


1 

1 

1 

L 

i 

i 

1 

1 
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1 

1 

I 
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] 

t 


0 D.A d.,+D.A 

A step by step breakdown of the various integrals can now be presented 
by simply shifting the filter envelope over the range of t". For clarity 
the amplitude of the filter envelope will be shown greatly enlarged. In 
addition, the time axis is not to scale. We define 

Z,(t) = a,D' cos(»„D't + J, ^ ^ 

There are four doppler dependent cases that must be examined in order 
to determine the exact values of the limits. 
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Case I. D' and D" are both less than 1. In this case the 


pulses are both expanded in time by the action of doppler. 

Case II , D' and D" are greater than 1. Pulse compression has 
occurred in this'instance. 

Case III . D' is less than 1, and D" is greater than 1. 

Case IV . D' is greater than 1, and D" is less than 1. 

Figure A-1 will be repeated and the relative positioning of the 
filter envelope will be shown to determine the proper limits in each 
instance. 


Case I. 


A. -T<t"<-(T-D^A) 



0 T+t" 

tJ^(T",W,A) 



(t) dt 



t 


B. -(T-D.A)<t"<0 
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C. Depending on the relative values of doppler, the next three integrals 
divide into different regions of t". 

For the case where D^A<^i-T 

Cl. 0 <t"<D^. 



Dl. D.A<t" J,-T 



El. I,-T<t"<I+D^-T 






































If, however, D.A is greater than ^i-T, then the following integrals 
result. 



D2. I,-T<t"<D.A 



E2. D^A<t"<J„+D^.A-T 



The remaining integrals in Case I do not depend on the particular 
value of D^. 
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T/D' T/D"+D.A 

/Z, (t) dt + /z (t) dt 

i" tti ^ 



i|<(T",a),A) 


T/D"+D.A 

/zjt) dt 

t" 


! 

j 
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The procedare is the same in the remaining cases. Consequently, 
only the integral expressions resulting from investigating these limits 
will now be presented; 

Case II. 


A. -T<t"<D.A-T 


B. D.A-T<t"<I-T 


ij»(T",a),A) = 


i|»(t“,w,A) = 


T+t" 

/z,(t) dt 

0 


T+t" T+t" 

/z.(t) dt + dt 

0 D. 


Dll + D^.A>T 


Cl. 


D1. 0<t"<J„+D.A-T 


V;(t",u,A) = 
i(<(t",w,A) = 


T/D' 


T+t" 


/z,(t) dt + /z^Ct) dt 

0 IJ^A 

T/D' T+t" 

J[z,(t) dt + dt 


D^A 


El. J„+D.A-T<t"<D.A ij;(T",a),A) = 


Dll + D^.A<T 


C2. J,-T<t" J>.+D.A-T i(>(t",w,A) = 


D2. ■^„+D.A-T<t"<0 


E2. 0<t"<D^.A 


ij;(T",w,A) = 


\j<(T",a),A) = 


T/D' T/D"+D^.A 

/Z,(t) dt + /Z,(t) dt 
t" CflA 


T/D' T+t" 

/z.(t) dt + dt 

0 D^A 

T/D' T/D"+D^A 

^Z,(t) dt + dt 

T/D' \/D"+D^A 

/z,(t) dt + /Z2(t) dt 

t' ILA 


F. D^A<t"<J, 




T|<(T",a),A) = 


tp(T",a),A) = 


T/D' T/D"+D.A 

/ Z,(t) dt + dt 

t" t" 

T/D"+D.A 

dt 
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